rrr=mmreader('C:\Users\georg\Documents\FFOutput\_MouseOnBallWithPS~1.avi');
vidFrames = read(rrr);
for k = 1 : numFrames
    mov(k).cdata = vidFrames(:,:,:,k);
    mov(k).colormap = [];
end

mmm=mov;

ps_ind=[51:75 174:199 262:287];


for ind=ps_ind
    mmm(ind).cdata(110:170,40:100,1)=255;
    mmm(ind).cdata(110:170,40:100,[2 3])=0;
end

movie(hf,mmm,1,25)



for ind=1:length(mmm)
    eee(:,:,ind)=mmm(ind).cdata(:,:,1);
end



figure;
plot(mean(ps_snips(:,ps_snip_time(:,10)<-0.01 & ps_snip_time(:,4)==4),2)-mean(mean(ps_snips(1:1000,ps_snip_time(:,10)<-0.01 & ps_snip_time(:,4)==4))),'g')
hold on
plot(mean(ps_snips(:,ps_snip_time(:,10)<-0.01 & ps_snip_time(:,4)==1),2)-mean(mean(ps_snips(1:1000,ps_snip_time(:,10)<-0.01 & ps_snip_time(:,4)==1))),'k')

axis tight
box off
xlim([500 6000])
set(gca,'xtick',[0:1000:6000])
set(gca,'xticklabel',[-1:5])
xlabel('Time [s]')
